Implicit ladder summation in the Hartree-Fock-Bogoliubov approach 
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The fully variational Hartree Fock Bogoliubov approach for bosons is studied in the limit of zero 
range forces in two- and three-dimensions. The equation of state obtained in two-dimensions is 
expressed in a parametric form. It is shown that the A potential permits to perform an implicit 
summation of the ladder diagrams without leaving the variational scheme, restoring thus the con- 
sistency of this approximation. 
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I. INTRODUCTION 

Modeling interatomic forces by a point-like isotropic 
potential appeared fruitful in the context of ultracold 
atoms as a way to study the dilute and low temperature 
regime for bosons or two-spin component fermions. In 
this respect, the Bethe Peierls model where pairwise in- 
teractions between particles are replaced by contact con- 
ditions is especially relevant to finding non trivial proper- 
ties in the two- and few-body problems, like confinement 
induced resonances (Tj-ll] or also in universal collisional 
properties as in the four-fermion problem Q. Neverthe- 
less the Bethe Peierls formalism is not adapted to usual 
formalisms of the many-body problem where the inter- 
action is described by a standard two-body vertex 0-Q]- 
The Fermi-Huang pseudo-potential which encapsulates 
the Bethe Peierls contact condition is a way to imple- 
ment directly the zero-range limit in a perturbative ap- 
proach of the many-body problem and was successfully 
used in the description of the dilute Bose gas [8|. How- 
ever when applied to the fully variational Hartree Fock 
Bogoliubov (HFB) approach, the Fermi-Huang pseudo- 
potential describes a molecular Bose gas instead of the 
usual atomic phase This problem is solved by using 
the A potential which in exact treatments, provides a de- 
scription strictly equivalent to the Bethe Peierls model, 
while it permits to improve approximate formalisms (Tol j . 
This potential is a function of the free parameter A which 
offers the opportunity to suppress the gap in the HFB de- 
scription of the homogeneous Bose gas, making thus the 
formalism compatible with the Hugenholtz Pines theo- 
rem [ll|. In Ref. [Hj], the formalism was applied to the 
two-dimensional (2D) Bose gas. It provides an equa- 
tion of state (EOS) which coincides with the result of 
Popov in the dilute limit [U, El] while giving at finite 
areal density an approximation which compares nicely 
with exact Monte-Carlo computations and perturbative 
results [l5l - fT7j ]. Moreover the formalism solves the prob- 
lem of the unexpected stabilization of a vortex in a non 
rotating trap, a phenomenon which occurs when HFB 
equations are truncated by neglecting the effect of the 
anomalous average of the noncondensate field operator 
[l9l [20l | . However, the derivation of the HFB equations 
using the A potential is not straightforward because the 



ansatz made for the density matrix doesn't belong to 
the correct Hilbert space (it doesn't support all the non- 
trivial singularities linked to the Bethe Peierls contact 
condition) . 

The present paper gives a simple derivation of the fully 
variational HFB equations for the three-dimensional (3D) 
and the 2D Bose ga s by using the A — e potential in- 
troduced in Ref. [21j . This potential enables to relax 
the strict zero range limit for the variational calculation 
while keeping the structure of the original A potential 
and the exact zero range limit is achieved afterward. This 
way, HFB equations are obtained without any approxi- 
mation following closely the standard and general varia- 
tional scheme described for example in Ref. [6[ . Emphasis 
is put on the 2D case where a parametric form of the non 
trivial EOS is derived. Finally, the ladder approximation 
of the Bethe Salpeter equation is solved in the low en- 
ergy limit. This calculation shows that the particular 
choice made for the parameter A in the HFB equations 
permits to implicitly sums the ladder diagram in the di- 
agonal self-energy at a collisional energy consistent with 
the evaluation of the off-diagonal self energy. 



II. INTERACTION POTENTIAL OF 
VANISHING RANGE 



A. Contact condition 



This paper deals with dilute gases of neutral bosons 
moving in a D-dimensional space (only the cases D = 2 
or D = 3 are considered) where interactions between par- 
ticles can be described by a short range pairwise potential 
of radius denoted by bo- At sufficiently low density (de- 
noted n) and temperature, i.e. rib D -C 1 and Xdb 3> 
where Xdb is the de Broglie wavelength, the transition 
matrix in two-body collisions of particles in the gas can 
be approximated by its lowest order in the binary colli- 
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sional energy (denoted E): 



TZa [ . ] is a regularizing operator: 
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where fco is the relative momentum (E = -r^f-), M is the 
reduced mass, au is the D-dimensional scattering length 
and 57 £) is the full Z5-dimensional space angle (f2 3 = At: 
and Q 2 = 27r). In Eq. ([T]) in the two-dimensional case, 7 
is the Euler's constant, and the 2D scattering length a 2 is 
by definition a positive parameter. The zero range poten- 
tial approach is a way to exactly recover these expressions 
of the transition matrix. It can be formulated for a sys- 
tem composed of iV-particles as follows. The many-body 
wave-function in the configuration space is singular at 
the contact (r^ = r 2 ) of each pair of interacting particles 
labeled by (ij) while out of these singularities it solves 
the Schrodinger's equation without pairwise interacting 
potential. For each interacting pair (ij) as the relative 
coordinates (denoted r^) tends to zero the singular be- 
havior is parameterized by the scattering length with: 
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+ 0( nj ); (D = 3) 



•0(r tf ); (D = 2). 



, (2) 

In Eq. @ the limit is taken for fixed values of (k 7^ i,j) 
and of the center of mass of the pair (ij) and the 
strength's singularity denoted by A L ^ 3 is a function of 
these coordinates. In this formulation often called Bcthc 
Peierls model in the literature, the pairwise potential is 
thus replaced by asymptotic conditions on the wavefunc- 
tion. While its simplicity is appealing, the Bethe Peierls 
model is not adapted to the standard techniques used in 
many-body theories. 



B. Pseudo-potential 



Another way to implement the zero-range potential ap- 
proach is to use the A-potential introduced in Ref. [lCj . 
In this formalism, the true finite range pairwise potential 
is replaced by a pseudo-potential which is in fact a way 
to perform implicitly a subtractive regularizing scheme 



achieved at the negative energy E A 



h A 

2ft, 



where A is 



an arbitrary positive parameter |21| . It can be introduced 
by its action on a two-body wavefunction ip as: 



(ti,t 2 \Va\iI>) =T D (E A )KAbP]6(r 12 ), 



(3) 



where Tjj(E a ) is evaluated at negative energy in Eq. ([T]) 
by using the usual analytical continuation ko = iA and 



11a 



lim 

ri2->0 



(d ri2 + A) [r n ■ } ; 
[l-ln(f Ar 12 ) r 12 d ri2 ] 



(D = 3) 



(D = 2). 
(4) 

In Eq. Q , the subtraction cannot be performed at A = 
in 2D while in 3D the A potential coincides with the 
Fermi-Huang pseudo-potential for A = 0. 

In Ref. [21|, it has been shown that in two- and few- 
body problems it can be fruitful to relax the strictly zero 
range limit by replacing the (S-distribution in Eq. @ with 
a family of functions 8 £ of finite support of range e while 
keeping the subtractive structure of the A potential and 
to perform the zero range limit e — > in a second step. It 
appears convenient to choose a family of Gaussian func- 
tions such that lim^oi 1 "!^) = S(r) which have the inter- 
esting property in the momentum representation to have 
an expression independent of the dimension D: 



(k\8 e ) = Xe(k) = cxp 



(5) 



Using this family of functions, the A — e potential can be 
written in the momentum representation as: 



(ki,k 2 l^ A lka,k4 



(2i:) n S(K 12 - K 34 ) 



T D (E A ) X Ak 12 ) limr A [ X£ (k 34 ) ■] . (6) 

e— >0 

In Eq. (J6)) k; is the momentum of the particle i, Ky 
(respectively k^) is the total (respectively the relative) 
momentum of the pair (ij) and is a regularizing op- 
erator: 



|A,«. 



(D = 3) 



(7) 



1 + - ln(e 7 A 2 e 2 /2)a e J - j ; (15 = 2). 



One may recognize in Eq. ([6]) a form similar to a separable 
potential. This structure is preserved in the configuration 
space: 



(r 1 ,r 2 |K A |r 3 ,r 4 )[-] = <5(R 12 - R34) 

T D (E A )(r 12 \5 t ) limr A [(«5 £ |r 34 ) 

£— >0 



(8) 



where R^- (respectively r^ ) denotes the center of mass 
(respectively the relative) coordinates of the pair (ij) and 
\r\8 e ) = (2ire 2 )- D / 2 exp (~r 2 /e 2 ) . 

It is worth pointing out that the way the zero range 
limit is performed in the A — e potential in Eqs. (|6|8j) 
totally differs from the usual method of the renormaliza- 
tion of the coupling constant (which tends to zero in the 
limit of a large ultraviolet momentum cutoff) often used 
in the many-body problem, even if the resulting equa- 
tions coincide when the zero-range limit is well defined. 
This issue is considered in the appendix 1X1 
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III. HFB FORMALISM FOR BOSONS 
A. Variational approach in the zero range limit 

The asymptotic conditions in Eq. ([2]) define the domain 
of the Hilbert space in the zero-range limit. In principle, 
in a variational approach one must preserve this structure 
by using an ansatz in the correct Hilbert space. However 
in the many-body problem, this constraint appears too 
stringent for being adopted. In this respect, the introduc- 
tion of the A — e potential greatly simplifies the imple- 
mentation of the variational schemes: one performs the 
variation of the functional at finite e avoiding thus all the 
difficulties linked to the singularities of Eq. and the 
zero range limit where e = is only achieved when one 
uses the equations resulting from the variational princi- 
ple. In the present section, this techniques is applied to 
the Hartree Fock Bogoliubov approach. 



B. Variation of the functional 

Using the expressions of the A — e potential in 
Eqs. (|6I8|) the HFB equations are obtained at finite e 
along the same lines as in Ref. |6j. This section gathers 
the major steps of this derivation. 

The gas is composed of N identical bosons of mass m. 
The field operator denoted by ?/> verifies the usual bosonic 
commutation rules: [?/>(r), -(/^(r')] = (2Tr) D 5(r — r') and 
[ip(r),i/j(r')] — 0. The standard HFB approximation is a 
symmetry breaking approach where the atomic field ij) is 
split into a classical field $ and a quantum fluctuation 



= %jj - $ with, (<j>) = 0. 



(9) 



The density operator is chosen as a Gaussian ansatz in 
terms of the field d>: 



exp(—K e /k B T) 



(10) 



where Z e is the partition function and the variational 
Hamiltonian K e is quadratic in terms of {(/),(/)'}: 

Ke = J d D ri d D r 2 [/i £ (r 1 ,r 2 )0 t (r 1 )^(r 2 ) 

+ iA e (r lj r a )0t( ri )£t( r2 ) 

+ iA:(r 1 ,r 2 )0(r 1 )0(r 2 )]. (11) 

This ansatz is an implicit function of the classical field 
$. The choice of a Gaussian density operator permits 
to use Wick's theorem. Consequently all the averages of 
products of the operator cj>(r) can be expressed in terms 
of the diagonal one-particle density matrix 



and of the off-diagonal one-particle density matrix 

/s e (ri,ra) = (0(r 2 )0(ri)) 6 . (13) 

In Eqs. (|10H13p and in what follows, quantities or av- 
erages where an index e is added means that they are 
considered for a finite value of e while the index e is 
dropped when the zero range limit (e — ► 0) is achieved. 
For example, this is the case for the depletion denoted 
by n(R) = /5(R, R) and for the total particle density 



n(R) = |$(R)| 2 +n(R). 



(14) 



In the subsequent equations the pair wavefunction is 
used: 

Ke (r 1) r 2 ) = (Vi(r 2 )V;(r 1 )> e (15) 
and also the A-regularized anomalous averages 

k a (R 12 ) = limr A / d^ru^ul^/s^rx.ra) (16) 

e^O J 

k a (R 12 ) = $(R 12 ) 2 + £ A (R 12 ). (17) 

In Eq. {jTTJ) k A can be also evaluated from Eq. (|4]) 
by h A = TZa [k] ■ The variational principle applied 
to the Grand-potential E e — TS e — p,N e where E e is 
the total energy, N e is the number of particles and 
S £ = — fcB(hi(i5 e )) e is the entropy gives @: 



& e (ri,r 2 ) = -7T- (A r <S) (ri 2 ) - jlS(r 12 ) + V C xt(Ri 2 ) 
+ 2 J d D r 3 d D r i (r 1 ,r 3 \V e A \r 2 ,r i }p e {r il r 3 ) (18) 

and 

A e (r l7 r 2 ) = J d D r 3 d D r i (r u r 2 \V e A \r 3 ,r 4 ) Ke (r 3 ,r i ). 

(19) 

After a straightforward integration one obtains: 



/i e (ri,r 2 ) = -— (AS) (r 12 ) + [V cxt (R 12 ) - fl] 6(r 12 ) 
2m 

+ SS n (R 12 )(r 12 |5 e ) (20) 



and 



A e (n,r 2 ) = ft£ 12 (R 12 )(r 12 |5 e 



(21) 



where the diagonal and off-diagonal self-energies are de- 
fined by: 



SE 11 (R) = 2T D (£ A )n(R), 
nS 12 (R) = T D (E A ) K A (R). 



(22) 



^e(ri,r 2 ) = (0 t (r 2 )0(ri)) £ 



(12) 



The equation verified by the condensate field is also ob- 
tained easily at finite e by minimization of the Grand 
potential. In presence of an external potential 14 x t(r) 
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acting on each particle of the gas, in the zero range limit 
it reads: 

{H] p - #$(r) + T D (E A ) [2n(r) + |$(r)| 2 ] $(r) 

+ T D (S A )^ A (r)$*(r) =0. (23) 

In Eq. ([23)) . the one-particle Hamiltonian H^ p is defined 
by: 



2m 



A r + V ext (r). 



(24) 



The variational ansatz K e is diagonalized by using a spe- 
cial Bogoliubov transformation on the particle field: 



5^ [6„u n (r) + SJX(r) 



(25) 



where the bosonic operator 6„ (resp. bjj destroys (resp. 
creates) a quasi-particle of quantum number n. The 
modal amplitudes {u n ,v n } in Eq. (|25|) are the eigenvec- 
tors of the linear system: 

(H lp - jj, + KSn) u n + hY,i 2 v n = hjj n u n 

(H lp - ft + Tffin) v* + HT,* 12 u* = -£u>„<. (26) 

In Eq. (|26p the coordinates (r) have been omitted for 
conciseness. In this basis the variational ansatz can be 
written as: 



K e = y^haJnbjbn. 



(27) 



From the commutation relations of the operators 
(b n , b' n , <j>, t ), one can deduce the closure relations 

K(ri)<(r 2 ) - <(nK(r 2 )] = <5(r 12 ) (28) 

n 

K(r!)<(r 2 ) - <(nK(r 2 )] = 0. (29) 

n 

The two fields p and k can be expressed in terms of the 
mode amplitudes at thermal equilibrium : 



p(ri,r 2 ) = J2 [<( r iH( r 2) (1 + fn) + 

Th 

u„(ri)<(r 2 )/„] 

K(ri,r 2 ) = ^y^ + fnj K(ri)un(r 2 ) 4 

<(r 2 )u„(ri)], 

where /„ = (b^ n b n ) is the Bose occupation factor: 

1 



fn 



«p(&*o-i' 



(30) 



(31) 



(32) 



In the zero range limit, the one-particle density p(r 1 ,r 2 ) 
is regular at ri 2 = and is equal to the depletion, while 



the pair wavefunction «(ri,r 2 ) verifies the Bethe-Peierls 
contact condition in Eq. @ as the relative coordinates 
r\2 tend to zero. The simplest way to show this prop- 
erty is to derive the equation verified by /«(ri, r 2 ). Using 
Eqs. (|26l28H3ip one obtains: 

[Hll + H\l\ K (r l7 r 2 ) + (r l7 r 2 |U A | K ) = 

2/2K(ri,r 2 ) - [fiSn(ri) + E u (r 2 )] K(r x ,r 2 ) 

- ^S 12 (r 1 )p(r 2 ,r 1 ) - ^S 12 (r 2 )p(r!, r 2 ). (33) 



In the right hand side of Eq. (|33l) there is no 6- 
distribution while the A potential acts on the pair wave- 
function in the left hand side of this equation. Hence, the 
pair wavefunction verifies the contact condition Eq. ([2]) 
and using Eq. ((4]) one founds the strength's singularity 
in terms of the anomalous average: 

2/i 



A K (R) = 



h 2 n 



-T d (E a )k a (R) va > 0. 



D 



(34) 



The off-diagonal self-energy evaluated from Eq. (|22|) is 
then proportional to A K : 



ftEia(R) 



2/t 



A K (R) 



(35) 



and is thus not an explicit function of the parameter A. 



C. A judicious choice for A(R) 

In an exact approach, the description of the system 
does not depend on a particular choice made for the 
free parameter A [which can be also a function of Ri 2 
in Eq. ©]. However, the expression of the diagonal 
self-energy in Eq. (|2"21 shows that HFB equations are A- 
dependent. Interestingly, this makes possible to solve the 
so-called gap problem in HFB (or Gaussian variational 
approach) which is not compatible with the Hugenholtz- 
Pines theorem [ll[ and was first shown to occur in 
Ref. [l8j |. For this purpose, the natural choice for the 
free-parameter A is given by imposing that the anoma- 
lous mode (uq, vq) = (<f>, — $) is a zero energy solution of 
Eq. (|2l))) . In an inhomogeneous systems, this leads to a 
value of A which depends on the position. This specific 
function A(R) is denoted by A*(R) and is obtained by 
solving the equation: 



k K ~ (R) = V R. 



(36) 



Using the A freedom of the off-diagonal self energy, one 
finds 



T D (£ A *)(R) =T D {E A )(R) 



1 



k (R) 
$ 2 (R) 



In the 3D case Eq. fl37|) gives 0: 

k a=0 (R) 



aA*(R) 



(37) 



(38) 



K A =°(R) +$ 2 (R)' 

while in the 2D case there doesn't exist such a general 



expression. 
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D. Homogeneous gas 



3D Bose gas 



In what follows equations are used for the homoge- 
neous gas with the choice A = A*, assuming thus that 
Eq. (|36[) is satisfied. The thermodynamic limit is also 
supposed to be achieved. The bosonic field operator is 
expanded on the plane wave basis: 



ip(r) = c + 



d D k 



2*0 



rCke 



(39) 



where $ = (co) and Ck annihilates a plane wave of mo- 
mentum k. The operators Ck, c^, verify the usual commu- 
tation rules: [ck, c k ,] = (27r)- D <5(k — k') and [c^, <v] = 0. 
The special Bogoliubov transformation takes the form: 

(2tt) 3 V 

where the operators (6k, b k ) respectively create and de- 
stroy a quasi-particle of momentum k and verify also the 
bosonic commutation rules [&k> b k ,] = (2Tr) D S(k — k'), 
[bk, &k'] = and the modal amplitudes verify 



& k it k e ikr + b k v k e 



-ik-r 



(40) 



\Vk\ 2 = l. 



(41) 



Equations (|23I36|) give the expression of the chemical po- 
tential for the homogeneous gas in terms of the conden- 
sate density and of the depletion: 



fi = T D (E A .)(<f> 2 + 2h). 



(42) 



Hence using Eqs. (|22|42|) one obtains for A = A* the re- 
lation p, = hT,n — ft£i2, which is the Hugenholtz Pines 
theorem [ll| . The solutions of Eq. (|26p are chosen real: 



Vk 



hbJu 



1 

V2 



h 2 k 2 

2m 



T D (E A *)& 



fVJJk 



1 

71 



K 2 k 2 

2m 



T D {E A *)<$> 2 



h 2 k 2 \* 



2m 



h 2 k 2 

2m 



- 1 



2T D (E M )& 



(43) 



At zero temperature the depletion is given by 

^\vk\ 2 (44) 

and in the zero-range limit, the off-diagonal one-particle 
density matrix is: 

f d D k —hYj\2 

K(ri ' r2 W (2^F^r exp(zk ' ri2) - (45) 

HFB solutions are found for a given value of the particle 
density: 

d D k „ 



n = 



(27T) 



D v k- 



(46) 



The 3D homogeneous case has been studied in 
Ref. PH. At zero temperature, Eqs. (|38|43I46[) together 
with the following expression of the regular part of the 
pairing field [obtained for example by using Eq. (|A1[) ]: 



~A=0 



12 



d 3 k 



(2tt) 3 \H 2 k 



1 



2Hojl 



lead to 



(n-$> 3 



30? 



[a 3 (3n- 2$ 2 )]' 



(47) 



(48) 



Hence, for each value of the atomic density, the conden- 
sate fraction is found by solving equation (|4"8j) . In the di- 
lute regime (no? < 1) one has n ~ <l> 2 , T 3 (E A *) ~ T 3 (0) 
and as expected the variational approach coincides with 
the usual perturbative Bogoliubov theory |l 0] . The tran- 
sition matrix T 3 (E A * ) is an increasing function of the 

density and diverges at the critical value n c = ^ 3 

where A* = 1/a and $ = 0. In this regime the mean 
field approach breaks-down, aA* is close to unity and 
the probability density for a quasi-particle to have a 
wavenumber k [i.e. v k k 2 / '(2 r K 2 n)\ is peaked near k = 1/a: 
collisions between non condensed atoms are not well de- 
scribed in this formalism. Clearly, inelastic processes and 
correlations not included in the HFB approximation must 
be taken into account to study this regime. 



2. 2D Bose gas 

In the 2D homogeneous case, the standard HFB for- 
malism for bosons can be used only at zero temperature 
where a true condensate exists. This situation was stud- 
ied in Ref. [12[ where it was shown that the present HFB 
equations permit to obtain an EOS which appears re- 
cently to be very good when compared to exact pertur- 
bative results beyond the Bogoliubov approximation. In 
what follows, this EOS is derived in a parametric form. 
For this purpose, one first calculates the asymptotic be- 
havior of the off-diagonal one-particle density matrix as 
r\2 tends to zero: 



K(ri ' r2) = ^T ln 



~yJmT 2 {Et,*)\<S>\ x ru 



which provides the A-regularized part: 

~ mT 2 (E A *)\$\ 
h 2 A 2 



K = A b U1 



0(r 12 ) 
(49) 

(50) 



The consistent choice for A* obtained from Eq. (|36p is 
thus 



A* 



/ -2tt|<I>| 2 
ln(eTA*a 2 /2)' 



(51) 
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The depletion is given by 

. _ mT 2 (_E A ,)|$| 2 



A 



* 2 



(52) 



4nh 2 An 
Finally using Eq. (|42l) , the constraint A = A* imposed by 
the Hugcnholtz Pines theorem provides the EOS 



t 2 [1-2 * ft*)] 



2,2 



4nra$ 



1 



(53) 



defined through the dimensionless parameter t = 02 A* 
(0 < t < 2e~ 7 ). The energy per particle is obtained from 
E/N = f ™ jjdn 



E 
N 



2 + 2 



hi 



v 2 ) 



1 



21n(^) 



(54) 



The low areal density regime corresponds to the limit of 
vanishing values of the dimensionless parameter t — > + 
and in this regime one recovers the result of Popov [l3| : 

4ft 2 



4irh 2 



In 



for na 2 "C 1, 



(55) 



e (27+l) m ^ a 2 

which is also the prediction of the Bogoliubov theory [3] ■ 
For finite areal densities, Eq. gives an estimate com- 
patible with Monte-Carlo simulations [H|,[l6j], the pertur- 
bative expansion one order beyond the Bogoliubov theory 
[l7l | and another approximate method (22f . 



where q = (qo,q). In what follows, this equation is solved 
at zero temperature by using the one-body Green's func- 
tion consistent with the HFB approximation: 



G(k) = 



(59) 



The structure of Eq. (1551) together with the definition of 
in Eqs. (|6I5T[) . give the form of the dressed vertex: 

(fci,fc 2 |r £ |fc 3 ,fc 4 >[-] = (2tt) D (5(K 12 - K34) 

3cffXe(ki2)limr £ A [ Xe (k 3 4)-] (60) 

and one wishes to compute the strength g e g. Integration 
over go gives: 



i f dq 



hi 2tt 



Gfe - q)G(k 2 + q) 



2 2 
W ki-q W k 2 +q 



2 2 
"ki-q^+q 



— Ml + ftw kl _ q + ftcj k2+q Ml + ftw kl _ q + fiWk,+q ' 

" (61) 

where — ui\ + In the limit of vanishing energy: 
k% —> , fc 2 — > and 51 — > one has: 



ffeff = 



l-T D (E A )J> 



(62) 



5. j4n implicit ladder summation 



where 



In terms of Feynman diagrams, the diagonal self- 
energy fi.Sn in Eq. (|2"2"|) is treated at the one-loop level 
[7J. That is why it explicitly depends on A [unlike the 
off-diagonal self-energy as shown in Eq. (|3"5j) ]. This sec- 
tion shows that the choice A = A* permits to implicitly 
perform the summation of the ladder diagrams enabling 
thus to evaluate the T-matrix appearing in SEu at the 
correct energy. 

In the standard many-body technique, the collision be- 
tween two quasi-particles is described by the effective 
potential T, which takes into accounts the effect of the 
medium Q: 

tei.&irife,,**), ( 56 ) 

where k = (w,k) denotes a D + 1 dimensional vector. 
In subsequent equations, the bare vertex for the A — e 
potential is denoted by: 



(k 1 ,k 2 \T^\k 3 ,k 4 ) 



(ki,k 2 |v; A |k3,k4 



(57) 



The ladder summation is obtained through the ladder 
approximation of the Bethe Salpeter equation Q: 



(ki>b.2\^e\k3i hi) — (kitk^lT^k^,^) + 



d D+1 



h J (2tt) d + 1 
x (fci,fe|r*|k-£,k + £)G(fci-£) 

x G{k 2 + q){k 1 -q,k 2 + q\T e \k 3 ,k 4 ), (58) 



J = — lim r. 



dD <l f N 2 "q + W q 

Xe{q) 



(27T) 



D ■ 



One then may use the identities: 



2u\v\ 



and UqWq 



(63) 



(64) 



The bit associated with M q v q in Eq. (|63|) gives an un- 
physical infrared divergence due to the fact that the 
effective potential is evaluated here in the strict zero 
energy limit of the scattering processes. In what fol- 
lows this contribution on the value of J is denoted 
by Jdiv A natural infrared cut-off for Jdw consistent 
with the ladder approximation in the medium is given 
by the inverse inter-particle distance q m i n ~ n 1 / . Ow- 
ing to the inequality <7 m ; n ^> l/£ where £ = y/ hj (2mSi2) 
is the healing length in the Bose gas one finds that 
T D (E A )J div ee C>(g-f n £- 6 ) can be neglected in Eq. 
The approximation it q + u q ~ 1 can thus be performed, 
and one may recognize in Eq. (|63|) the expression of the 
A-regularized part of the off-diagonal one-particle density 
matrix at zero temperature considered for finite values of 
r: 



K e (ri,r 2 ) 



(27T) 



2hu!r, 



■ exp(zq • n 2 ). (65) 
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One thus obtains J = k A / (HT,^), and using Eqs. (|37I62[) 
the strength appearing in the effective potential is given 
by: 

g eff = T D (E A ) $2 = T D (E M ). (66) 

Assuming that the spectrum is gapless, this result is valid 
for all A. For A = A* then T e = 1^ in Eq. flSSJl, meaning 
that when used at the one-loop order the bare vertex 
incorporates the ladder summation for A = A*, which is 
the desired result. 



where the coupling constant is adjusted to reproduce the 
low energy scattering properties is well adapted for nu- 
merical studies and leads to equations very similar to the 
one obtained using the contact conditions of Eq. [23[ . 
This appendix focuses on a standard techniques of the 
renormalization of the vertex interaction in the 3D con- 
tinuous space where the strength of the force vanishes 
in the zero range limit. It is shown below that in the 
zero-range limit of the effective interaction, one recov- 
ers in the two-body problem the contact condition of the 
Bethe Peierls model obtained from the A — e potential in 
the momentum representation. 



IV. CONCLUSIONS 

In this paper the variational HFB equations in the zero 
range potential approach have been obtained following 
the standard derivation of Ref. [6]. This way, results 
of Refs. [HI EH where the A potential was used are re- 
covered in a simple manner. The EOS for the 2D Bosc 
gas obtained in the HFB approximation is expressed in a 
parametric form. It is explicitly shown that the specific 
choice made for the parameter A permits to perform the 
ladder summation while keeping the variational scheme 
intact. These results have been obtained by using a rep- 
resentation of the A potential where the strict zero range 
limit is relaxed in intermediary calculations. The natu- 
ral continuation of this work is to use the A — e potential 
vertex of Eq. ([57]) in standard techniques of many-body 
theory borrowed from quantum field theory, including 
the treatment of the fermionic gas. The A freedom is a 
gift of the zero-range potential approach which permits to 
improve approximate treatments of the many-body prob- 
lem. The best choice for the parameter A in the general 
time dependent HFB formalism or also in the variational 
Balian Veneroni approach [6J remains an open question. 
Another interesting issue is the possible use of a fully 
variational approach for the 2D Bose gas at finite tem- 
perature where the low energy classical modes of the field 
operator could be treated consistently in the variational 
scheme. 
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1. Bethe Peierls contact condition in the 
momentum representation 

In Ref. [2l|, the operator r A was introduced without 
loss of generality through the identity: 



lim r' 

£^0 



lim 

e-s-0 



h 2 n 



D 



2//. 



(Al) 



where \ip 6 ) is a two-body state considered in the center of 
mass frame, is its strength's singularity in the limit 
where e tends to zero, and \4>^) is a reference state equal 
to the action of the two-body Green's function at the 
energy E A = — ^77" on the ket \S e ): 



2fj, 



A 2 ' 



(A2) 



Equation (IA1|) illustrates the fact that the A potential is a 
regularization techniques based on a subtractive scheme. 

Using the operator , the Bethe Peierls contact con- 
dition in Eq. @ can be written without referring to any 
representation as: 



lim r* = 



n D h 2 



A. ir 



2fi T D (E A 



(A3) 



Inserting the closure relation in the left hand side of 
Eq. (|A3|) . the contact condition can be expressed in the 
momentum representation as: 



lim 

e-s-0 



d°k 



(27T) 



D 



Xe(k) 



(k|* e 



k 2 + A 2 



/u(iA) 
(A4) 

where A > 0, /d(«A) = -2^T D (E A ) / '(^d^ 2 ) is the scat- 
tering amplitude in the Bethe Peierls model evaluated at 
negative energy and 



Am = 



-1 

£ln k 



lim (fc 2 lim(k|# £ 

:— s-oo V e — >0 



(A5) 



Appendix A: Contact condition obtained from the 
renormalization of the coupling constant 

In the literature, numerous effective interactions are 
used for implementing the zero-range potential limit. As 
an example, the contact interaction on a discrete lattice 



For example, in the A — e potential model a scattering 
state of incoming momentum k reads: 



(k|Ve 



where A. 



(27r) 3 <5(k - k ) 
-fn(k ). 



k'o 



k' 2 



(A6) 
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2. Link with the renormalization of the coupling 
constant 



In what follows, Eq. (jA4J) is derived by using the model 
potential: 



V(r) = v Q g Q {r), 



(A7) 



where vq is the bare coupling constant and ffQ(r) is a 
family of functions converging toward a ^-distribution in 
the limit of arbitrarily high values of the ultraviolet mo- 
mentum cutoff Q (Q — > oo). For a two-body scattering 
state with an incoming plane wave of momentum |k ), 
the transition matrix of the model is 



(k|T(i; + *o + )|ko) = (k|v|* ko > 



(A8) 



The Fourier transform of the function gg(r) is denoted 
by 



5Q(k) = J d 3 rexp(-ik • r)g Q {r) 



(A9) 



so that the transition matrix reads: 



(k|r(£ + iO + )|ko) =v Q 



d 3 k' _ 



(2tt) 3 

A simple choice for the function gQ is 



fa(k-k')<k'|*k >. 

(A10) 



9Q 



1 if k < Q 
otherwise 



(All) 



In the limit of small momenta as compared to the ul- 
traviolet momentum cutoff Q (|k| <C Q and |ko| -C Q) or 
equivalently in the limit of large ultraviolet momentum 
cutoff, the following identity holds: 



d 3 k' 

— <?Q(k-k')(k'|* k0 ) 



d 3 k' 

Q^°°Jk'<Q ( 27r ) 3 



k'|* k0 ), 



(A12) 

which makes possible to directly identify the scatter- 
ing amplitude of the model at low collisional energy in 
Eq. (TATOj) : 



/q(ko) 



2irh 2 
Q— >oo 2-rth 2 



k |r(£^o + )|k ) 



k'<C 



d 3 k' 

(2^ 



(k'ltfi 



(A13) 



Comparing with Eq. (|A6I) one obtains the relation 

d 3 k' 



Aq, = lim 



q-s-oo 27rft 2 Jk'<Q (2?r 



(k'|*k„ 



(A14) 



Using the fact that 



<k|^) = (2^(k-ko) + <k 5^t <0 ty , (A15) 



E + i0+ - 

one obtains a closed equation for /<g(ko): 



— / (ko) = l-47r/ Q (ko) 



2/i 



d 3 fc' 



fc'< 



q (2tt) 3 Jfeg + i + - fc 2 ■ 



(A16) 



The integral in the right hand side of this last equation 
is performed easily by using the standard analytical con- 
tinuation fco = iq (where q > 0) and finally: 



/q(ko) 



1 



27th 2 j_ 2Q 

HVQ 



^+ik 



(A17) 



The value of the bare coupling constant vq is chosen such 
that the scattering amplitude /g(ko) coincides with the 
Bethe Peierls scattering amplitude /s(fco) in the limit of 
an arbitrarily large momentum cutoff Q, thus yielding: 



1 

«3 



2nh 2 



2Q 

IT 



(A18) 



Using the fact that | a.3 1 <C Q, the bare coupling constant 
can thus be expressed as: 



7T 2 h 2 



7T 3 H 2 



V Q 



From Eqs. (|A18[) and (|A14|) one obtains: 



(A19) 



lim 



k<Q 



d 3 k 

(2^)3 



(k|*t 



k 2 



a3 



(A20) 



Finally, the identity 



lim , 

Q^JkKQ (27r) 3 



d 3 k ( 47T 47T 

k 2 k 2 + A 2 



= A (A > 0) 



(A21) 

permits to rewrite Eq. (|A20|) in the same form as the 
contact condition of Eq. (|A4|) . 
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